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Abstract 

The rate per unit volume for anomalous electroweak baryon number violation 
at high temperatures, in the symmetric phase, has been estimated in the 
literature to be 0(a^T 4 ) based on simple scaling arguments. We argue that 
damping effects in the plasma suppress the rate by an extra power of a w to 
give 0(a^T 4 ). We show how to understand this effect in a variety of ways 
ranging from an effective description of the long-distance modes responsible 
for baryon number violation, to a microscopic picture of the short-distance 
modes responsible for damping. In particular, we resolve an old controversy 
as to whether damping effects are relevant. Finally, we argue that similar 
damping effects should occur in numerical simulations of the rate in classical 
thermal field theory on a spatial lattice, and we point out a potential problem 
with simulations in the literature that have not found such an effect. 
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I. INTRODUCTION 



Anomalous baryon number (B) violation in the hot, symmetric phaseQof electroweak the- 
ory occurs through the creation of non-perturbative, nearly static, magnetic configurations 
with spatial extent of order ((jf 2 T) _1 , where g is the electroweak SU(2) coupling. (This is 
reviewed below.) It has been widely assumed that g 2 T must then be the only scale relevant 
to the problem, so that the baryon number violation rate per unit volume, by dimensional 
analysis, must be 0((g 2 T) 4 ). We shall argue that damping effects in the plasma cause the 
time scale for transitions through these configurations to be 0(g 4 T) instead of 0(g 2 T), and 
therefore the rate is 0((g 2 T) 3 (g 4 T)) = 0(a^T 4 ). 

The possible importance of damping effects was pointed out many years ago in ref. |§ 
in the context of B violation in the symmetry-broken phase of the theory. The effect was 
controversial, and another analysis of the problem |4j] claimed that damping plays no role. 
In the intervening years, many people have privately expressed the concern that damping 
might affect the symmetric phase rate, and we make no claim to be the only, or even the 
first, people to think of it. The purpose of this paper is simply to resolve the controversy, 
elucidate the physics involved, and put to paper the result that the rate is 0(a^T 4 ). Many 
of the individual parts of this paper will be reviews of various aspects of thermal physics 
already familiar to some readers, but we believe their synthesis in this discussion of the 
baryon number violation rate is original. 

In the remainder of this introduction we briefly review the conventional picture of anoma- 
lous transitions in the symmetric phase and the standard estimate of the rate. Then we give 
a quick but formal estimate of the effects of damping that closely follows the broken-phase 
discussion of damping in ref. [|J. 



1 Here and throughout, we use the term "symmetric phase" loosely since, depending on the details of 
the Higgs sector, there many not be any sharp boundary between the symmetric and "symmetry-broken" 
phases of the theory . By symmetric phase we shall mean temperatures high enough that the magnetic 
correlation length is 0(\/g 2 T) and determined by non-perturbative dynamics. 
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The body of the paper is devoted to providing a number of different ways in which to 
understand the physics underlying this result. Section II begins with an analysis of the 
power spectrum of gauge field fluctuations based on the fluctuation-dissipation theorem 
and argues that it is fluctuations with frequency of order g A T and spatial extent (^ 2 T) _1 
which have sufficiently large amplitude to generate non-perturbative transitions. This is 
followed by a real-time, finite-temperature diagrammatic analysis which demonstrates that 
perturbation theory breaks down in exactly this domain due to the non-linear interactions 
of low-frequency, low-momentum components of the non-Abelian magnetic field. Then, in 
section III we sketch the microphysical origin of damping and turn to one of the traditional 
methods of understanding the transition rate: by computing the rate at which the system 
crosses the "ridge" separating classically inequivalent vacua. We show that a single net 
transition from the neighborhood of one valley to the next actually involves 0(l/g 2 ) back- 
and-forth crossings of this ridge, and so standard methods which count individual ridge- 
crossings overcount the true rate by a factor of 1/g 2 — precisely the suppression due to 
damping. In the final section, we discuss how damping effects should also play a roll in the 
topological transition rate for classical thermal field theory on a spatial lattice. This appears 
inconsistent with the results of numerical simulations 0, which show no sign of damping 
effects. We explain why the simulations done so far may fail to measure the true topological 
transition rate. 

This multiplicity of approaches in not intended to be a case of several poor arguments 
"substituting" for one good one. We believe that any one of the viewpoints discussed below 
provides compelling evidence that thermal damping is responsible for suppressing the high- 
temperature electroweak baryon violation rate by one power of a w relative to the naive order 
ct%T^ estimate. Our goal is simply to examine the phenomena from multiple perspectives 
and understand the connections between them. 
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FIG. 1. A schematic representation of the (bosonic) potential energy along a particular direction (labeled 
z) in field space, corresponding to topologically non-trivial transitions between vacua. 

A. Review of the standard picture 

Fig. |l] is the standard visual aid for thinking about anomalous transitions. Consider 
the theory in Hamiltonian formalism or in A = gauge, where the degrees of freedom are 
A(x) and the conjugate momenta are E(x) = —dA(x)/dt. The horizontal axis represents 
one particular direction in the infinite-dimensional space of gauge configurations A(x). The 
minima represent the vacuum A = and large gauge transformations of it, labeled by 
their Chern-Simons number N cs . The vertical axis represents the potential energy of the 
configurations.^ Whenever a transition is made from the neighborhood of one minimum 
to another (which we call a topological transition), the electroweak anomaly causes baryon 
number to be violated by an amount proportional to AiV cs :Q 



2 More precisely, the non-fermionic contribution to the potential energy. When a transition is made, there 
will also be the perturbative energy cost of the fermions created by that transition. 

3 There are also other directions in configuration space along which N cs changes that have nothing to do 
with the vacuum structure of the theory and exist even in U(l) theories. These directions are not ultimately 
relevant to baryon number violation. To make the issue more precise, imagine starting with a cold system 
with some baryon number, heating it up for a time to make anomalous transitions possible, and then quickly 
cooling it. The system will cool into the nearest vacuum state shown in fig. [l]. So the net change in baryon 
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That said, we shall now ignore the fermions and, for simplicity, focus on the rate for topo- 
logical transitions in the pure bosonic theory. And, since we are interested in the symmetric 
phase of the theory, we shall generally ignore the Higgs sector as well and focus on pure 
SU(2) gauge theory.^ 

Topological transitions can occur at a significant rate if thermal fluctuations have enough 
energy to get over the top of the energy barrier separating neighboring vacua. The most 
important parameter determining the energy of the barrier is the spatial extent R of the 
configurations depicted in fig. |l|. In order to generate an 0(1) change in Chern-Simons 
number (or baryon number), Eq. ( |1 . 1| ) plus dimensional analysis implies that a gauge field 
of spatial extent R must have a field strength^ of order (gR 2 )^ 1 . Consequently, the energy 
of these configurations on the potential energy barrier will be 0((g 2 R)~ 1 ), while the gauge 
field amplitude itself (smoothed over the scale R) is 0((gR)~ r ). 

So a more representative picture of the configuration space is provided by fig. 0a: the 
energy barrier is not a single point but a ridge which becomes arbitrarily low if arbitrarily 
large configurations are used to cross it. The smallest R for which non-perturbative thermal 
transitions through such configurations are not significantly Boltzmann-suppressed is R ~ 
l/g 2 T. For the same reason, this is also the scale where perturbation theory breaks down 
in the hot plasma. 



number in this example depends on whether the system has made a transition from the neighborhood of 
one vacuum state to the next, not simply on whether Chern-Simon number has temporarily changed due to 
an excursion in some irrelevant direction. See Q for a more detailed discussion. 

4 Except in the immediate vicinity of the electroweak phase transition (or crossover), fermions or Higgs 
fields merely provide additional "hard" thermal excitations in the symmetric phase and do not affect any of 
the following discussion in a substantive fashion. 

5 We are using the conventional gauge field normalization in which no factor of 1/g 2 multiplies the kinetic 
terms in the action and perturbative fluctuations are 0(1) in amplitude. 
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FIG. 2. The same as fig. [l] but supplemented by an extra dimension of configuration space corresponding 
to (a) spatial size R of the configurations, and (b) some other generic direction y, such as a particular mode 
of momentum ~ T . 

R ~ l/g 2 T is in fact the dominant spatial scale for topological transitions. One argument 
is entropy: there are fewer ways to cross the barrier with larger configurations than smaller 
ones. Another argument is due to magnetic confinement in the hot plasma. Ignore E = —A 
at the energy barrier for the moment, so that the configuration is purely magnetic. Though 
static non-Abelian electric forces are no longer confining at high temperature, static magnetic 
forces are.0 The confinement scale is just the spatial scale of non-perturbative physics, l/g 2 T. 

Now consider the transition rate. Suppose that a configuration of size R on the energy 
barrier ridge of fig. ||a is produced in the plasma. What would be the time scale for its 
decay? The unstable modes of the configuration will be associated with momenta of order 
1/R. If the configuration were decaying at zero temperature, the time scale for decay would 



6 The standard way to see this is to consider the expectation of very large spatial Wilson loops at high 
temperature. These can be evaluated in the Euclidean formulation of finite-temperature gauge theory, 
where Euclidean time has a very small period /? at high temperature. By dimensional reduction this is 
equivalent to understanding the behavior of large Wilson loops in three-dimensional Euclidean theory. But 
three-dimensional SU(2) gauge theory is confining. 
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then be 0(1/ R). So the standard estimate is that there is one unsuppressed transition per 
volume R 3 per time R, giving the rate T ~ R~ A ~ (g 2 T) 4 . As we shall see, the flaw in 
this estimate is that the configuration is not decaying at zero temperature but instead in 
interacting with the other excited modes of the plasma. 

We should emphasize that the picture of fig. ^|a is still incomplete because we have 
shown only two degrees of freedom in configuration space. There are an infinite number of 
other degrees of freedom in which the potential turns up so that, for fixed R, the energy 
barrier looks like a saddle as depicted in fig. |]b. The energy ridges discussed before are 
really hypersurfaces through configuration space, separating valleys that contain the different 
vacua. One can formally define these hypersurfaces as follows 0]: For each configuration 
A(x), follow the steepest descent path away from it in the infinite-dimensional analog of 
fig. |2| — that is, follow (minus) the gradient of the potential energy. For generic configurations, 
one will eventually approach one of the classical vacua. The barrier ridge hypersurface which 
separates the classical vacua is the exceptional hypersurface which doesn't flow to a vacuum 
configuration. We shall return to this picture of the ridges in detail later, when we take up 
the old controversy about whether damping can affect the transition rate. 

B. A quick estimate of damping 

Consider again the decay of a configuration Ab(x) of size l/g 2 T on the energy barrier. 
As the system passes through Ab, one may analyze the motion by linearizing the equation 
of motion. Expanding A(x, t) = Ab(x) + <5A(x, t), the equations have the form 

d 2 6A{x) = J(x) - K 2 5A(x) + 0(5 A 2 ) , (1.2) 

where J is the gradient of the potential along the ridge and K 2 is the potential energy 
curvature operator with a negative eigenvalue — k, 2 of order (g 2 T) 2 corresponding to decay 
away from the ridge. For the moment, let's focus on that one eigenmode, reducing our 
equation to 



S 



d 2 SA = k 2 5A + 0{5A 2 ). (1.3) 

The solution to this equation grows exponentially on a time scale of l/g 2 T. 

Now consider the interaction of this decay with typical thermal excitations of the plasma. 
The transition rate involves physics at soft energies g 2 T that are small compared to the hard 
energies T of typical particles in the system. The simplest way to analyze the problem is 
to consider an effective theory for the soft modes of the theory, where the physics of the 
hard modes has been integrated out. In the context of (|1.3|) , we need to know the effective 
(5 A) 2 interaction generated by integrating out the hard modes. The soft-mode self-energy 
generated by the hard modes is dominated by the processes shown in fig. || and is known as 
the hard thermal loop approximation to the self-energy. 7 



Q hard Q hard Q hard Q hard 



^ ^ 

Psoft Psoft Psoft Psoft 

FIG. 3. The soft-mode self-energy generated by forward scattering off of hard modes. The external hard 
lines are on-shell. 



There is a reasonably simple formula for the resulting self-energy U(uj,p) in the limit 

u,p<CT ||, but for our purposes it will only be important to know its qualitative behavior 

in terms of the ratio rj = u/p: 

( g 2 T 2 (0(l)-tO( V )), 77 «1 , 
II L = (1.4a) 
g 2 T 2 0(l), 77 >1, 



7 The hard thermal loop approximation also generates corrections to cubic and higher couplings of multiple 
soft particles ||. These will be needed in a complete quantitative analysis, but they do not affect any of 
our simple estimates. 
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g 2 T 2 {0{rf)-iO{ri)), rj < 1 , 
n T = { (1.4b) 
g 2 T 2 0(l), V >1, 

where I1l and Ht are the longitudinal and transverse parts of the (retarded) one-loop self- 
energy. To distinguish them, it is convenient to momentarily switch to covariant gauge| 
(rather than A = gauge), where the longitudinal and transverse projection operators are 

~ -— /> - — ^ O - —) (1 5a) 

L ~ IpI 2 v p2 J v p2 J ' 
IpI 

where P' 1 = (u, p) and has only spatial components. The 0(g 2 T 2 ) behavior of I1l(0, p) 
reflects the Debye screening of static electric fields at distances of 1/gT. The vanishing of 
Ht(0,p) reflects the absence of similar screening for static magnetic fields. The 0(g 2 T 2 ) 
behavior of both IIt(u;, 0) and Hl(uj, 0) reflects the O(gT) mass gap for propagating plasma 



waves. 



Q hard Q-P hard 
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FIG. 4. Absorption or emission of soft modes by hard modes when to < p, contributing to ImU. 

The imaginary part of ELr arises from absorption or emission due to scattering off hard 
particles in the thermal bath, as depicted in fig. ||. By kinematics, these processes occur 
only when u < p. Now modify eq. (|1.3|) to include the self-energy (and Fourier transform to 
frequency space), 



8 For Ao = gauge, 7\ of ( 1.5a ) is restricted to spatial indices and is then —Pq/P 2 times the projection 
operator = p 1 ^ /\p\ 2 - There is a compensating factor of — P 2 /Pq in the longitudinal piece of the free 
propagator in A = gauge: i£> (0) = p- 2 (V T - (P 2 /Pq)Vl,). 
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-u 2 5A = (k 2 -U(tu))5A. (1.6) 

Given that k ~ p ~ g 2 T, the self-energy II (a;) will stabilize the unstable mode unless we 
focus on transverse modes and take 77 <C 1. We may then approximate 

n ~ ig 2 T 2 r] ~ zu;T , (1.7) 

so that the solution to the linearized Eq. (|1.6|) is 

uj~ig A T. (1.8) 

The decay time is therefore of order l/g A T, or 1/g 2 slower than assumed in the standard 
non-dissipative estimate. 

II. SOME ALTERNATIVE VIEWS 

A. Analysis in terms of spectral density 

We now consider another formal way to see that uj ~ g 4 T is the appropriate frequency 
scale for topological transitions. Start by considering the power spectrum pa(u, p) of gauge- 
field fluctuations A(u, p) in the plasma: 

(AV, p)*AV, P')> = pS(". P) (^) 4 5(lu-u') 5 3 (p- P ') . (2.1) 
This may be related to the retarded propagator 

iD l £(t-t',x-x!) = ([A^x^A^x')]) 6{t-t') 

= Z- l Tie-^ H {A\t^)\A j {t'^)}e{t-t') (2.2) 

by the fluctuation-dissipation theorem:^ 

9 The theorem follows from inserting a complete set of energy eigenstates in (^^), Fourier transforming, 
and then taking the imaginary part. See, for example, sec. 31 of ref. |l0[ . 
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Pa(uj, p) = -2(n w + 1) IhlDr(w, p) = -2(n w + 1) Im 



1 



, (2.3) 



(u + it) 2 - p 2 - U(u, p) 

where fi u = l/(e^ w — 1) is the Bose distribution function. 

Now, using this relation, consider which frequencies of A will have enough power to 
generate topological transitions. As reviewed in section 1, this requires fields with spatial 
extent R ~ ^/g 2 T and amplitude A ~ 1/gR ~ pT. Consider the right-hand side of 
and use the hard thermal loop approximation (|1.4j ) for II in 

ImIT 



lmD R (uj, p) 



(2.4) 



(u 2 - p 2 - Ren) 2 + (imn) 2 

The last term corresponds to propagating plasma waves. For p ~ g 2 T, the integrated 
power given by the right-hand side of (|2.3|) from all frequencies of order u is of order 

'l/g 2 T, tu-gT; 
uop A ~ ll/w, g A T<u<p ■ (2.5) 

. u/g 8 T 2 , u < g'T . 

For uj <C gT, it is dominated by transverse fluctuations. A schematic plot of the power 
Pa{w) is shown in fig. |||. From (|2.1| ), the power in a particular soft mode with R ~ l/g 2 T, 
and frequency of order uj, will produce fluctuations with amplitude 

' 9 2 T , u ~ gT ; 

g 3 ^/^, g A T<iu<p ; (2.6) 

g' 



A ~ (p 3 ujp A ) 



1/2 



VST, cj<^ 4 T 



Hence, non-perturbative amplitudes A ~ 1/ gi? ~ gT are generated when u; ~ g 4 T. 



10 It is only a ^-function in the leading-order approximation to II. The smearing of the 5-function by the 
plasmon width, and other features of the spectrum such as the two-plasmon cut, will not be important for 
the order-of-magnitude power estimates we make below. 
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FIG. 5. A schematic plot of the power Pa(w) of gauge-field fluctuations with p ~ g 2 T. The 5-function 
spike at lo ~ gT corresponds to propagating plasmons. 

In section 3, we shall see how to get this same power spectrum by considering the 
microphysical details of the behavior of the hard degrees of freedom, and we will translate 
the power spectrum into a qualitative discussion of what the actual real-time trajectories of 
the soft degrees of freedom look like. 



A pictorial way to represent the previous argument is to ask when diagrammatic per- 
turbation theory breaks down in the effective theory of the soft modes, since we need non- 
perturbative effects for a topological transition. For simplicity, let us focus on transverse 
modes with u <C p, for which one can ignore Reil. Consider adding a loop to a Feynman 
diagram, as shown in fig. ^ and consider the thermal contribution to that loop, shown in 
fig. ^|c. The cost of adding the loop is order (gp) 2 from the new vertices, p~ 4 from the two 
new uncut propagators, and d 3 pdu p^(u;,p) for the phase space probability of finding the 
new soft particle in fig. |6|c, where pA ~ n u Iul-Dr. Integrating over frequencies of order uj 
and momenta of order p ~ g 2 T, the cost of adding a loop is of order 



B. Estimate from Feynman diagrams 




g^T/u, g 4 T <u <p , 
u/g 4 T, uj < g A T . 



(2.7) 



So the loop expansion parameter is 0(1) for u ~ g 
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(a) 



(b) 



(c) 



FIG. 6. Adding a loop (a— >b) to a Feynman diagram, (c) represents the thermal contribution to the 
new loop, corresponding to forward scattering off a particle in the thermal bath. 



Instead of considering Feynman diagrams in the effective theory (in which the self-energy 
has been resummed), we can get some insight into the origin of the transition time scale by 
recasting the above argument in terms of diagrams in the original, microscopic theory The 
key to the estimates above was the behavior of ImilT, which arises from the interactions 
shown in fig. |j. The origin of the dominant contribution arising from multiple interactions of 
the form of fig. ^|c is interactions such as those shown in fig. |7|. The straight lines represent 
hard thermal particles; the wavy lines represent virtual magnetic quanta with (u,p) ~ 
(g 4 T, g 2 T) that are absorbed or emitted. 



C. Interacting magnetic fields 




FIG. 7. A microphysical picture of the important interactions for topological transitions. 
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The cost to the rate of adding a new interaction with a new hard particle is (p~ 4 ) 2 for 
new propagators, (gp) 2 for a new soft vertex, (gT) 2 for a new hard vertex, 0(1) for hard 
particle Bose or Fermi factors, and 

P-P^- 0(|q-q'| S p) 0(K-ovl £ w) ~ (2-8) 
for phase space. The total cost, by this estimate, is 

(p- 4 ) 2 x (gp) 2 x (gT) 2 x 0(1) x Tup 2 ~ (2.9) 

and is indeed 0(1) for ~ g 4 T. For larger lu, where this microscopic perturbation theory 
appears to go completely wild, one must consider the effect of resumming the effects of Imll 
into the soft propagators and return to our previous estimates. The advantage to having 
considered fig. [7| is simply that it gives us a physical picture for the origin of topological 
transitions in the plasma. The soft, virtual lines emitted from the hard particles simply 
correspond to the low-frequency and low-momentum components of the magnetic fields 
produced by the movement of those particles. These are not propagating electromagnetic 
waves but simply the magnetic fields carried by all moving charged particles. Topological 
transitions are then created by the non-linear interactions of these fields. 



III. MICROPHYSICAL PICTURE OF BARRIER CROSSING 

The normal approach used to calculate the topological transition rate in the symmetry- 
broken phase, and to estimate it in the symmetric phase, is based on calculating the rate, 
in equilibrium, at which the system crosses the energy barrier hypersurface separating the 
classical vacua |3],f|]. If one assumes that each such crossing is associated with a net transition 
of the system from the neighborhood of one vacuum to the neighborhood of another, then 
the barrier crossing rate is a good measure of the topological transition rate. Implicitly 
making this assumption, ref. ||] claimed to show that damping has no effect on topological 
transition rates, in contradiction to the claims of ref. ||. In this section, we resolve this 
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dilemma by showing that the assumed equality of the barrier crossing and net topological 
transition rates fails if one studies the full, short- distance theory. The difference in the 
rates will turn out to be precisely the suppression uj / ' g 2 T ~ g 2 we have argued arises from 
damping. 

A. Review of the microscopic origin of damping and the Langevin equation 

Since the purpose of this paper is pedagogy and not excruciating details of formalism, 
we shall make some simplifications in order to elucidate the physics involved. First, we 
shall treat all the modes as classical and simply assume there is an ultraviolet cut-off at 
momenta of order T (which is where quantum mechanics enters to cut off the ultraviolet 
catastrophe of classical thermal statistical mechanics). Secondly, though we shall consider 
all the hard modes of the theory, we shall only focus on the one soft mode a(x) that is 
responsible for the decay of any particular configuration on the barrier. So we shall restrict 
6A(t, x) = A(t, x) - A b (x) to 

8A(t, x) ~ 5A soh a(x) + / d 3 q 8A(t, q) e iq x , (3.1) 

J hard 

where the amplitude of a(x) is normalized to 0(1) near its center. Finally, imagine putting 
the system in a box so that we can discretize the degrees of freedom. Rather than studying 
gauge theory directly, we shall begin by reviewing the derivation of damping for soft modes 
in a generic theory with cubic interactions: 

(5H = \fi + \K 2 z 2 + £ (\ Py f + n 2 \y t \ 2 ) +z^y*Gi j y J ■ (3.2) 

i ij 

Here z represents the soft mode, the hard modes, K 2 ~ —(g 2 T) 2 the curvature of the 
potential energy for the unstable soft mode [but we could also consider stable soft modes 
with K 2 ~ +(g 2 T) 2 }, Qi the hard frequencies of order T, and Q the soft-hard-hard part of 
the three-vector coupling. 

The basic approximation is to realize that, provided the coupling Q is perturbative, the 
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motion of the many hard modes is not much affected by the motion of the soft mode. To 
first approximation, then,Q 

y^a t e~^\ (3.3) 
where en, are random phases and thermally-distributed random amplitudes: 

(a*aj) = /(fii) 5ij = ^ 8 ij > ( 3 - 4 ) 



which is the equipartition theorem with our normalizations ( |3.2| ). [With continuum normal- 
izations, the quantity corresponding to f(fli) would be the classical limit T/fij of the Bose 
distribution n(Qi).} Now consider the equation of motion for the mode of interest, z: 

z + K 2 z = -y*Qy. (3.5) 

The leading approximation ( |3.3| ) produces an effectively random force term £(t) on the 
right-hand side, 

m = -E^ a ^~ m3 ' Q ' )t - ( 3 - 6 ) 

Using ( |3.4j ) , the time correlation of this force is0 

(£(t)£(t+At)> = \Gn\ 2 e-^-W* . (3.7) 

ij 

The force £(t) is therefore a source of colored noise for the evolution of the soft mode, and 



the power spectrum of that noise is given by the Fourier transform of ( |3.7| ) 



(£H£(u/)) = 2n6(u-u') ]T /(^)/(%) |^| 2 2^(^-0*-^) . (3.8) 



11 The simple time evolution ( |3.3| ) is to be understood as approximating the hard modes on time scales 
loss than their thcrmalization time. 
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Rigorously, the microcanonical problem involves once randomly choosing the on and then evolving the 



system, whereas ( 3A ) implies averaging over an ensemble of choices for each . As long as the number of 
hard degrees of freedom that interact with the soft mode is large, there is no essential difference between 
these two cases. 
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The square matrix element just represents the process of fig. |] and, if one returns to contin- 
uum language, the power spectrum turns out to be simply 

(£(u;, p)£(u/, p')> = -Att5(co-uj') <5 3 (p-p') n(w) ImllH . (3.9) 

So far we have reviewed the origin of an effectively random force term in the equation 
of motion for soft modes (a term that we did not discuss in our quick analysis of section 1). 
To see the origin of damping, one must consider the perturbative effect of the soft mode on 
the motion of the hard modes. The equations of motion for the t/i are 

y i + n 2 i y i = -z{t)J2G ij y j . (3.10) 



The perturbation to the solution is easily obtained by using the free solution ( |3.3j ) for y on 
the right-hand side and then Fourier transforming: 

yi(uj) = ai2ir6(uj-tti) + - — — ^- — V^w-fi^i/yc^ + 0(z 2 ) , (3.11) 

(uj + ley — iq j 

where we have used the retarded solution for the response of yi to x. Putting this solution 
into the soft-mode equation ( |3.5| ) now gives 

(-co 2 + K 2 )z = + 2 J2 GijGj* a*a k ±^hz^l + (f) . (3.12) 

ijk 10 1 h 

Averaging the second term on the right-hand side over the random amplitudes ati finally 
yields 

{—lo 2 + K 2 + H(uj))z ~ , (3.13) 

where 

nM=2p (n ,) (M + n| l %,_^ . (3.14) 

This is just the discretized version of the soft-mode self-energy of fig. [| 

The purpose of this review has been to show that the effective equation of motion ( |1.6| ), 
which we first used to argue that the decay time is 0(l/g 4 T), should more accurately be a 
Langevin equation of the form 



[-uj 2 + K 2 + 1%)) 6A soit = £(u) , (3.15) 



where £(u;) is a random force with (£) = and power spectrum ( |3.9|) . For a stable soft mode 
(K 2 > 0), the decay of the mode due to the imaginary part of the U(u)5A term and the 
excitation of the mode by the random force term balance each other to maintain thermal 
equilibrium of the soft mode. 

Now (dropping the subscript "soft" ) consider the solution 

SA(t) = 5A(t) + AA(t) , (3.16) 

where 5A(t) is the solution to the homogeneous equation that we originally considered in 
section 1 and AA(t) are the fluctuations induced by £: 

{ -^ + V + m a) y (3 - 17) 



Computing the power spectrum of AA using (3.9|), one recovers our earlier result of (|2.3| 



for the power pa, projected onto the soft mode under consideration. 



B. What does the decay look like? 



Let's now return to our characterization (|2.5|) of the power of the fluctuations in A 
integrated over frequencies of order uj. We'd now like to investigate how much the relevant 
soft mode oscillation AA(t) wiggles in time so that we can assess whether the barrier is 
crossed once or multiple times per net transition. 

First, how many times does the motion of AA(t) change direction per unit time? This is 

equivalent to asking about the fluctuations in AA(t). We can obtain the integrated power 

in A simply by multiplying (|2.5|) by uj 2 : 

T , uj ~ gT 5 

up A ~ < u , g A T <uj < p ; (3.18) 

, u 3 /g 8 T 2 , uj < g'T . 

Unlike the power for the amplitude A of fluctuations, the power for A is dominated by plasma 
oscillations u ~ gT and not by low frequencies uj ~ g 4 T. The time scale for the motion of 
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AA(t) to change direction is therefore 1/gT. The power in AA and AA can be summarized 
by considering three characteristic frequency scales: gT, g 2 T, and g A T. The results of ( f2.5|) 
and ( |3.18 ) are summarized for these scales in table | and depicted schematically in fig. |8|. 
Remember that the oscillations AA are superimposed on top of the slow net motion of the 
homogeneous solution A for the damped decay. 



OJ 


\A\ 


\A\ 


gT 


g 2 T 


g 3 T 2 


g 2 T 


g 2 T 


g'T 2 


g 4 T 


gT 


g 5 T* 



TABLE I. Amplitude of fluctuations in A(t) and A(t) corresponding to three characteristic frequency 
scales. 



AA(t) 



EH 




l/g 4 T 



FIG. 8. A schematic picture of the time evolution of the fluctuations AA(t) showing fluctuations of scales 
oj r~j gT, g 2 T , and g T. Keep in mind, however, that there is really a spectrum of fluctuations between 
oj ~ g 2 T and oj ~ g A T whose amplitude grows as uj decreases. 



Now let us estimate how many times the system crosses 5 A = during a topological 
transition. The time scale for the net transition is t net ~ l/g 4 T, over which 5 A changes 
magnitude by S net A ~ gT. During that time, 5A(t) changes directions of order gTt net ~ l/g 3 
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times due to plasma oscillations whose amplitude is 0(g 2 T). The spectrum of fluctuations of 
the system may be rather complicated, but the total time the system spends within 0(g 2 T) 
of S A = will be of order 

*netX-^-~— -. 3.19 

b nct A g 4 T 

During this time, the plasma fluctuations are large enough to drive the system back and 
forth across 8 A = with frequency gT. Thus, during one net transition, the system crosses 
through 5 A = of order 1/g 2 times. Any calculation which simply computes the rate of 
barrier crossing per unit time will overestimate the rate of net transitions by a factor of 

i/g 2 - 



C. Being more precise about the "barrier" 

In the previous section, we implied that crossing SA so f t = many times during a transi- 
tion implies that the energy barrier hypersurface is crossed many times during a transition. 
However, 5A so f t is just the projection of A — Ab onto the soft mode of interest, and the 
exact equation for the energy barrier hypersurface is not really as simple as 5A so { t = 0. 
Perhaps, as the hard modes oscillate, the 5A so f t location of the barrier also oscillates. In 
order to confirm our picture, we should check that the effect of hard modes of the shape of 
the barrier is small enough not to affect our argument. We verify this in Appendix [A]. 



IV. IMPLICATIONS FOR CLASSICAL LATTICE THERMAL FIELD THEORY 

The estimates of damping we have made have all been based on continuum, quantum, 
thermal field theory, where the effective ultraviolet cut-off for any thermal effects is 0(T), 
since the distribution of particles with momenta ^> T is Boltzmann suppressed. One of the 
numerical testing grounds used in the literature for studying topological transition rates, 
however, has been classical thermal field theory on a spatial lattice. In classical thermal field 
theory, modes of arbitrarily high momenta are thermally excited (leading to the ultraviolet 
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catastrophe of the classical blackbody problem), and the ultraviolet cut-off is the inverse 
lattice spacing a -1 rather than T. As we shall argue below, this infinite growth in the 
number of relevant hard modes as a — > leads to infinitely strong damping of the transition 
rate in the continuum limit of classical thermal field theory. 

In the limit u, p <C T, the usual continuum result for the imaginary part of the transverse 
self-energy Ht from figs. [3] and f| is of the form 

IhiITt(cc>, p) ~ g 2 J ^A n g - \M- W \ 2 ] <*Hq| - p ■ q) 



2 w 



~ 9 - / dqqn q 
p Jo 

-g 2 T 2 -, (4.1) 
p 

where M.^ is the vertex in the left-hand figure of fig. [| The q integral is dominated by the 
ultraviolet and cut-off by q ~ T for the quantum case, shown above. 
For classical field theory, the only difference is that 

n Q -> — (4.2) 

and the ultraviolet cut-off is now the inverse lattice spacing a -1 , so that the last step of ( |4.1|) 
now gives 

Imn c i(a;,p) ~ ^Ta" 1 -, (u < p) . (4.3) 

p 

Damping will therefore suppress the transition rate by an extra factor of aT compared to 
the quantum field theory case, so that the rate is 0(a 5 aT 5 ) and vanishes in the a — > limit. 

This result is in apparent contradiction with the numerical results of ref. H, which 
claimed to find 0(a 4 T 4 ) behavior for the rate. Our result is suppressed by an additional 
factor of 0(aaT), which is one power of the dimensionless lattice coupling (3. Ref. f| used 
values of j3 ranging from 10 to 14, so they should find a 40% violation of the assumed 0(a 4 T 4 ) 
scaling of the rate. Such an effect is clearly inconsistent with their statistical errors. (See 
their fig. 3.) 
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There is a possible problem, however, with the assumption that ref. || is in fact mea- 
suring the topological transition rate. The problem stems from the subtleties of trying to 
define topological winding number on the lattice. The idea in ref. || was to measure a 
lattice analog of the square of the topological winding: 

(J*dt J d 3 xtrFF^ = ((AN cs (t)f) , (4.4) 

where it is now convenient to consider field strengths normalized so that the action is F^ u /g 2 . 
Under the picture of fig. [I], the topological winding in pure gauge theory should randomly 
diffuse away from zero, so that at large times (fOj) behaves like TVt where T is the transition 
rate. The authors implement this procedure on the lattice by making a lattice approximation 
to trFF, which schematically has the form of the cross-product 

trFF -> E a x ti(Ua a ) = J , (4.5) 

where U represents plaquettes perpendicular to the links the electric field E lives on. (See 
ref. || for the detailed expression.) The potential problem with this representation is that, 
unlike the continuum expression for trFF, the integral (i.e., lattice sum) of J over space is 
not a total time derivative. The time integral of this lattice analog to / cPxtrFF does not just 
depend on the initial and final configurations but depends on the path used to get from one 
to another. In particular, consider a path in time that starts from some configuration near 
the vacuum, never makes any non-perturbative excursions from it (and so in particular never 
makes a topological transition), and finally ends up back at the initial configuration. The 
lattice analog of the left-hand side of fl4.4|) , (/ dt d 3 x J) 2 , would not be zero. Perturbative 
fluctuations can therefore either increase or decrease / dt d 3 xJ without increasing the energy, 
and so there will be a purely perturbative contribution to the diffusion.^ 

To estimate the size of this lattice artifact in the diffusion rate T, consider expanding 
( |4.5[) in powers of the lattice spacing a in lattice perturbation theory, remembering that 



13 This point has also been observed by Ambjorn and Krasnitz and will be addressed in a forthcoming 
publication. 
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U r^i e m 3 . The leading E ■ B term is a total time derivative and does not cause problems. 
As an example of a subleading term, consider a term involving E and three powers of B: 

trFF trFF + a^EBBB + ■■■ (4.6) 

Now consider the contribution of this term in the right-hand side of ( ]4.4|) : 

a 8 J* dh jT* dt 2 ^ (| d 3 x EBBB^j (^J d 3 x EBBB^j \ (4.7) 

The correlation is dominated by the ultraviolet and will be quasi-local in space and time. 
So, in the large time limit, it becomes 

a 8 Vt J dt' d 3 x ((EBBB) 0fi (EBBB) t , ^ ~ /T ~ 5 8 TVi (4.8) 

This lattice artifact therefore gives a contribution to the measured diffusion rate T of 
0(a T 4 ). The moral is that purely perturbative effects, having nothing to do with true 
topological transitions, might obscure the true topological transition rate, which we have 
argued is 0(a 5 aT 5 ). 

This work was supported by the U.S. Department of Energy, grant DE-FG06-91ER40614. 
We would like to thank Jan Ambj0rn, Patrick Huet, Alex Krasnitz, Larry McLerran, Misha 
Shaposhnikov, and Frank Wilczek for useful conversations (some of which date back ten 
years) . 

APPENDIX A: HARD MODE EFFECTS ON THE BARRIER SURFACE SHAPE 

As discussed in section 1, the barrier hypersurface is the surface which (1) separates 
the vacua and (2) maps into itself when the gradient of the potential is followed. In this 
appendix, we want to focus on whether the hard modes have a significant impact on the 
shape of the surface. Specifically, consider (a) one soft, unstable mode of interest, and (b) 
all hard modes, as in the generic model of ( |3.2| ) with K 2 = —k 2 < 0. What is the equation 
of the barrier hypersurface in this subspace? 
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We shall restrict attention to typical hard mode amplitudes in the thermal bath, i.e., 
\yi\ ~ 1/fij. Secondly, we shall assume that the interaction Q between the soft mode and any 
individual hard mode is perturbative, as it indeed is hot gauge theory. The translation of 
this condition to the generic model ( |3.2| ) is easily made by considering stable soft modes, in 
which case \z\ is typically 1/k and our perturbative condition is xyQy <C 1, giving Q <C nVl? . 

Under these conditions, our result for the ridge equation is 

* = £tfn^fc = y^y- ( A1 ) 

This is easily checked as follows. First, it contains (z, y, y*) = (0,0,0) as it should. Next, 
we must check that 

W(z, y*, y) = (-k 2 z + y*£y, tl 2 y + zQy, y*tl + zy*G) (A2) 

lies in the tangent plane to the surface (within our approximations). On the surface (AT), 



W(z, y*, y) ~ (y^y, Q 2 y, y*&) . (A3) 
The tangent plane to the surface is spanned by 

-^(z(y*,y),y*,y) = ((Sy) h e h O) , (A4) 

^Wy*,y),y*,y) = ((^y)„o,e J ) , (A5) 
dyi 

and so 

n = (1, -Sy, -y*S) (A6) 



is normal to the surface. ( |A3|) and (A^) then give n* ■ W ~ as desired. 



To count the number of barrier crossings in section 3, we should really have studied the 
evolution of z — y*Sy rather than the evolution of z. We are now in a position to check 
whether this makes any difference by checking the amplitude of the fluctuations in y*Sy. 
Using the leading-order behavior ( |3.3| ) for the yi and averaging over the amplitudes yields 
the equal time amplitude 
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v (y^y) 2 ) = E/(^)/(^)l^l 2 - (A7) 

Comparing this to 

{e) = Ef(^)m 3 )\G l j\ 2 - (as) 



ij 

shows that ( |A7| ) is of order (£ 2 ) times the typical size of f2~ 4 (which is T -4 ). The continuum 
version is that the equation for the barrier surface is 

5 A = x(A hard (t)) (A9) 

where 

(X) = 0, (x 2 )~T- 4 (a (A10) 

and £ here is understood to be projected onto the soft mode. Using ( |3.9|) and p ~ g 2 T, 

(X 2 ) ~ T~V J duo n w Imn(cu) ~ T~Y x £ 2 T 3 ~ (g 4 T) 2 . (All) 

This 0{g i T) magnitude of the fluctuations in the location of the surface is much smaller 
than the magnitude g 2 T plasma oscillations described in section 3. The shape of the surface 
therefore has no effect on our previous discussion of crossing the barrier. 
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